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ABSTRACT 

Aims. We investigate the cosmological test based on the evolution of the gas fraction in X-ray galaxy clusters and the stability of the cosmo- 
logical parameters derived from it. 

Methods. Using a sample of distant clusters observed by XMM-Newton and Chandra, and comparing their gas fraction at different radii to the 
gas fraction observed for nearby clusters, we have determined the likelihood functions for £1„ in a flat universe and the confidence contours in 
the Q.,,, - Ha plane. 

Results. Results obtained at the virial radius point to a high matter density Universe, while for inner radii the Q.,„ parameter obtained tends to 
decrease, reaching values compatible with the concordance model. The analysis allows us to conclude that this test provides ambiguous results 
due to the complex structure of the ICM that induces a dependence of the gas fraction on temperature, radius, and redshift, which cannot be 
accounted for by the self-similar picture expected from pure gravitational heating of the ICM. 

Conclusions. The use of gas fraction in X-ray clusters to constrain cosmological parameters seems therefore to be compromised until a better 
understanding of the ICM physics and the ability to obtain observations of better quality up to the virial radius are achieved. 
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1. Introduction 

Physical proprieties of galaxy clusters have been considered a 
promising tool for cosmological probes. A test based on the as- 
sumption of a constant baryon fraction in these objects was pro- 
posed by |Sasaki (1996)| and |Pen (1997)| These authors showed 
that under this assumption, the gas fraction of clusters can, in 
principle, be used as a standard "ruler" to determine cosmo- 
logical parameters. The main idea is that if the only heating 
source of the gas is gravitational collapse, the baryon frac- 
tion should be constant and independent of redshift. Since the 
measured baryon fractions depend on the cosmology through 
the angular distance to the clusters, the "true" cosmological 
model would be one where all the measured baryon fractions 
at different redshifts would be equal. This method has been ap- 
plied in the past years (Rine s et al. 1999ilEtFori & Fabian 20001 
lAllen et al. 20021 . Ettori et al. 20031 lAllen et al. 2004 ) using 
clusters observed by Chandra and ROSAT PSPC, with the re- 
sults pointing to a low value of Q,„, consistent with the con- 
cordance model {Q.,„ x 0.3, Ha ~ 0.7). However, these authors 
have often restricted the method to a single radius, defined by 
a given high overdensity relative to the critical density of the 
Universe. To check the validity of this cosmological test, it is 
desirable to perform a good study on the behavior of the gas 
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fraction distribution (the main contribution of baryons) inside 
clusters and its evolution with redshift. Such a study was done 
recently in |Sadat et al. (2005)| hereafter SOS, who claim that 
clusters are more complex objects than it was expected in self- 
similar models, as they found that gas fraction presents a de- 
pendence on radius, temperature, and redshift. 

The goal of the present work is to give a quantitative analy- 
sis on the results obtained in SOS, by determining the likelihood 
functions of the cosmological test based on the gas fraction for 
the sample of SOS and study the effect of using radii of different 
overdensities to perform it. The paper is organized as follows. 
In Sect. 2, we present the cluster samples that was used in the 
test. In Sect. 3, we present the method used to determine the gas 
fractions for all clusters in the sample. The cosmological results 
are presented in Sect. 4 and finally, in Sect. S, we draw our con- 
clusions. We used a Hubble constant of SO km/s/Mpc, unless 
the dependence is explicitly given (with i/o=100/;km/s/Mpc). 

2. The cluster sample 

We used three samples of clusters taken from different refer- 
ences. The nearby sample consisted of 3S nearby clusters with 
redshifts between 0.02 and 0.09S, observed by ROSAT PSPC 
and compiled in |VikhIinin et al. (1999) This sample contains 
symmetric clusters with high-quality imaging data up to large 
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radii, which allows for good estimations of the gas masses 
on the outer parts of the clusters. We have also used XMM 
data obtained on a sample of eight high-redshift clusters ob- 
served as part of the XMM-Newton Q project, a systematic 
XMM - Newton guaranteed time follow-up of the most dis- 
tant SHARC clusters [(Bartlett et al. 200 1)[ The high sensitivity 
of XMM-Newton allows for an investigation of the emissivity 
in high-redshift clusters beyond half the viral radius, a remark- 
able result [(Arnaud et al. 2002)| The sample contains clusters 
within a restricted range of redshifts, between 0.45 and 0.65, 
with the detailed data reduction and analysis being presented 
in |Lumb et al. (2004)| 

To have a bigger statistical set of high-redshift clusters, 
we considered another sample of eight clusters observed by 
CHANDRA, within the same range of redshifts (0.451 to 0.583) 
as the above XMM clusters. This sample was taken from 
IVikhUnin et al. (2002)| who claim to have an accurate measure- 
ment of the X-ray brightness nearly up to the virial radius. 

For all clusters, we had information about the temperature, 
surface brightness parameter fi and core radius rc taken within 
the observed region of the clusters, and some information (de- 
pending on the sample) that allowed us to normalize the gas 
mass profile, all the ingredients needed to determine the gas 
mass fraction as it will be shown in the next section. 



3. Gas mass fractions 

3.1. Determination of the gas mass 

The gas mass fraction of a cluster at a given radius r is de- 
fined as fgasir) - Mgas{r)IM,oi{r), needing good estimations 
of the gas mass and the total mass contained inside a radius 
r. The gas mass is generally determined by integration of the 
so-called y6-model for the description of the X-ray gas profile 
|(Cavaliere& Fusco-Femiano 1976)[ 



Pgasir) = Pgaj(O) 



rc 



(1) 



which translates to the observed X-ray surface brightness: 



1+ - 
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The observed brightness is due to the well-known process of 
bremshtralung radiation, and so the central density p(0) can be 
determined from the central brightness 5'x(0), allowing a de- 
termination of the gas mass at any radius of detection with a 
good accuracy. It is important to notice that the obtained re- 
sults depend on the assumed cosmology because the conver- 
sion of the observed angles to physical radii depends on the 
angular distance to the cluster as r = 9dA{z\ O,,,, ^^a)- Using the 
y6-model we were then able to determine the gas mass profiles 
of our sample of clusters within the range of detection, for any 
given set of cosmological parameters. When necessary, the gas 
mass beyond the radius of detection and up to the virial radius 
was estimated by extrapolation using a y6-model with parame- 
ters obtained within the region of detection. 



3.2. Determination of tine total mass 

In the simplest picture of purely gravitationally driven forma- 
tion of virialized systems such as galaxy clusters, it is ex- 
pected that such objects exhibit self-similarity (Kaiser 1986), 
meaning that the radial profile of any physical quantity, such 
as density, should exhibit a similar shape that is independent 
of the cluster's reference mass and redshift, when normalized 
to the virial radius. Indeed, numerical simulations in which 
only gravitational physics is taken into account indicate that 
halos of diff'erent masses follow a universal density profile 
[(Navarro, Frenk & White 1997 )| the so-called NFW density 
profile. We considered this profile with a concentration param- 
eter of 5: 



p{r) 



= 1500 



'200 



r(5r + r2oo)^ 



(c = 5) . 



(3) 



The normalization of this profile can be obtained using the to- 
tal mass at the virial radius, which can be determined by its 
definition: 



Mv = -npoil +zy AvRl 



(4) 



In this expression Ay is the virial contrast overdensity com- 
pared to the mean density of the universe at redshift z and po 
is the mean density of the Universe today. Since many authors 
use a contrast overdensity Ac compared to the critical density 
of the Universe (pjz) = 3H(zf/S7TG), we found it useful to 
make our analysis using both definitions to compare with dif- 
ferent published results and confirm that results are indepen- 
dent of the adopted definition. With this second definition the 
last expression becomes: 



4 
3' 



Mv = -7Tpc(z)AcRl 



(5) 



The virial overdensity is usually computed from the spherical 
top-hat model for gravitational collapse, and for an Einstein-de 
Sitter universe, one has Ay = Ac - l^n^. We used the formulae 
given in |Bryan & Norman (1998)| which provide the method to 
determine Ay and A^. for any flat cosmological model. The scal- 
ing of the mass-temperature relation can be obtained from the 
expressions above, using the fact that for a purely gravitational 
heating of the gas, T oc GM/Ry, which yields: 



T = AtmM- 



2/3 



178 



1/3 



h^'\l+z)keW 



or. 



T = BtmM 



2/3 



{h'AcE^) 



1/3 



keV 



(6) 



(7) 



where M15 is the mass virial mass expressed in units of 15Mq 
and H(z) = HjE^ = i/2[Q„(l + z? + Qr(1 + zf + Qa]. 
The normalization constants Ajm or Bjm can be determined 
from numerical simulations ( |Evrard, Metzler & Navarro 1996[ 
Bryan & Norman 1998) and observations. Several studies have 
been performed to test the My - T relation as predicted from 
numerical simulations by means of X-ray observations (see. 
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Fig. 1. Likelihood functions of the gas mass fraction test at different reference radii, for a flat universe. Left panels: Assuming 



a constant normalizations Atm of 6.24 keV, (Bryan & Norman (1998) top) and 11.4 keV (down). Top right: Assuming a nor- 
malization Atm that depends on Q,„ as Atm = 4.9^"" '^ KeV. Lower right: Likelihood functions using normalization A7-m=1 1.4 
keV, but performing the test at the same physical radii as if Atm - 6.24. Solid lines correspond to the test applied at the virial 
radius. The other lines correspond to likelihood functions obtained for different raddi: Ri (dot), R2 (dashed), Rsoo (dashed-dot), 
and R2500 (dashed-dot-dot). 



e.g., Horner et al. 1999; Nevalainen et al. 2000; Finoguenov 
et al. 2001; Sanderson et al.2003; Pedersen & Dahle 2006). 
Some disagreement has been found concerning the normal- 
ization constant Atm and the slope (steeper than the pre- 
dicted 1.5) for cooler systems (Tx less than 4 keV). Recently, 
[Blanchard & Douspis ( 2005) considered an alternative way to 
determine Atm- By assuming that the observed baryon content 
of clusters is representative of the Universe's baryon content 
and using the WMAP tight constraint on Q./, and the relation 
between and h, they derived the mass-temperature normal- 
ization as a function of Q.,„, obtaining Atm - 4.9Q^,"'^^ + 10%, 
for 0.3 < Q.,„ <1.0. 

4. Cosmological results 

As it was mentioned, the use of the gas fraction of clusters as a 
cosmological test has been based on the fact that this physical 
quantity should be constant at a given radius when normalized 
to the virial radius. However, S05 claims that the baryon frac- 
tion of nearby and distant clusters presents a clear temperature 
dependence, which is larger in the inner parts and appears to be 
very shght at the virial radius. This fact indicates that, in prin- 



ciple, only the virial radius, if any, can be used as the reference 
radius to perform the cosmological test based on a constant 
gas fraction. One can relax this assumption, and consider that 
the dependence of the gas fraction on temperature for high red- 
shift clusters should be related through scaling arguments to the 
same dependence observed for low-redshift clusters. Indeed, if 
the gas mass fraction varies with T, scaling would imply: 

fg (R/Rv, T, z) = fg (R/Rv, T x r.(0)/r.(z), z = 0) , (8) 

where T,(z) is a characteristic temperature associated with a 
characteristic mass scale at epoch z and defined by cr(M»,z) a; 
5th- So, to use the observed gas fractions for local clusters and 
compare them to the distant ones, we proceeded as follows: for 
any given set of cosmological parameters, gas fractions were 
determined for the local sample of clusters, and a fit was made 
to determine an analytical expression for its temperature de- 
pendence in the form of fg{Rx) - Cx T"". This was done for 
five different normalized radii - the virial radius, two fiducial 
radii Ri and R2 (defined in |Vikhlinin et al. (1999)| as Riooo and 
R2000, with the index referring to contrast overdensities rel- 
ative to the background baryon density), and also two radii 
commonly used to analyze the gas fraction that correspond 
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to more central regions of the clusters - R500 and R25oo- For 
these radii, the index are referred to overdensities related to 
the critical density of the Universe. The results of the fits to- 
gether with Eq. dS) provide a model that allows the use of high- 
redshift gas fractions to determine the cosmological parame- 
ters, since the "true" cosmological model should allow us to 
have fg'^''"^(Rx, z, T) = Cx (r.(0)/r.(z))"^ T"*. The UkeUhood 
function for this test can then be written as ^-e\p{-x^/2), with 
being: 



X 



f /,,,(□,„, Qa, Zi) - fT'^Ti, zc, Q„, Qa) 



(9) 



In this expression, is the number of distant clusters in our 
sample and cTdisp is the dispersion around the best-fit of the lo- 
cal cluster distribution. The errors on the derived baryon frac- 
tion for individual clusters (cr,) were neglected since they are 
much smaller than (Jdisp- 

We used three different approaches to the normalization 
of the mass-temperature relation: two constant values Ajm - 
6.24 keV (Bryan & Norman 1998) and A™ = 11.4 keV, a 
Q„,-dependant value Ajm - 4.9Q,7*'^'^ (Blanchard & Douspis 
2005), and finally, we assumed Ajm - 11-4 keV, but taking the 
gas fractions at the same physical radii as if Aju - 6.24. This 
last approach was done in order to cross-check the dependence 
of the results on both the assumed radii (physical) to perform 
the test and the normalization constant. For the case of an Ajm 
dependent on Q„,, we did not make the analysis at R500 and 
R2500 because at these radii Eq. (7) cannot be applied since we 
do not have an expression for Bjui^m)- 

The obtained results are presented in Fig. 1. At first glance, 
the main feature in each plot is the fact that the likelihood func- 
tions are highly dependent on the considered radius: Q„, ob- 
tained from extreme radii are mutually exclusive and cannot 
therefore be combined to obtain a convergent value for Q„,. 
It is important to examine whether this result depends on the 
calibration of the mass-temperature relation, and we have ad- 
dressed this question in some detail: in the same figure we can 
see the effect of using different approaches to the normalization 
constant. After increasing Ajm by a factor of 2 (left panels), 
no significant changes are observed since the maximum likeli- 
hood values decrease only slightly. The same is observed when 
we used the Q„,-dependant normalization, where the results are 
similar to the case of Ajm - 6.24 keV. This is consistent with 
the fact that changing the normalization modifies mass in a 
similar way for low and high redshift clusters. Finally, when 
we considered a value of 11 .4 keV for Ajm, but taking the gas 
fractions at the same physical radii as with Ajm - 6.24 ke V, the 
likelihood functions remain almost unchanged, except for the 
virial radius. At this specific radius we can see a great decrease 
in the maximum likelihood value, since the radius is now much 
larger than the actual virial radius. This shows that this test is 
less dependent on the adopted M - T normalization at inner 
radii, although even at this radii we cannot obtain compatible 
likelihoods for Q„,. 

We can point out that, as expected, no difference was found 
in the results at the radii studied by 805 (Ri,,>, Ri, and R2) 
when using the two alternative ways to determine the virial ra- 
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Fig. 2. The confidence contours for Q,„ and Qa at the 1, 2, and 3 
cr level on 1 parameter using the cosmological test based on the 
gas fraction of distant and nearby clusters and assuming Ajm - 
6.24 keV. Top: Test performed at R2. Bottom: Test performed at 
the virial radius. 

dius and mass (comparing to the mean or the critical density 
of the Universe). We have also determined the confidence con- 
tours in the two-dimensional parameter space of Q„, - Qa. For 



this, we used the formulae provided by |Pierpaoli et al. (2001) 
to calculate the virial overdensity Ay for any set of cosmo- 
logical parameters and followed the same analysis described 
above. Results are presented in Fig. 2 for the test at the virial 
radius and at R2. The contours are compatible with the likeli- 
hood functions obtained in the case of a flat Universe and they 
clearly show the significant difference in the results of the gas 
fraction test when we consider an inner or outer radius to per- 
form it. 

5. Discussion 

We have used combined data from Chandra and XMM-Newton 
to apply the cosmological test based on the gas fraction evo- 
lution considering that this quantity varies with the tempera- 
ture of the cluster and that this dependence should evolve with 
redshift through a simple scaling argument based on the gravi- 
tational growth of structure in an expanding Universe. The ob- 
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tained results allow us to conclude that this kind of test is highly 
dependent on the radius considered to implement it, since we 
found significantly different values for when the test is per- 
formed at different normalized radii. The complexities in clus- 
ter's internal structures which were evidenced in 805 are at the 
origin of this observed ambiguity. It seems that the gas frac- 
tion observed in low and high-redshift clusters cannot be ac- 
counted for by a purely gravitational heating of the gas, which 
in other words means that non-gravitational effects present in 
clusters induce a significant deviation of the gas fraction pro- 
file from standard scaling. We found no significant differences 
in the results when using two different definitions for normal- 
izing the radius, and it is interesting to point out that the results 
obtained at Rasoo are compatible with those obtained by Allen 
et al.(2002) and Allen et al.(2004) at the same radius. Note al- 
though that these authors did not apply the same method, as 
they assumed a constant gas fraction at their working radius. 

The main question raised by this work is whether there 
is a radius at which scaling relations can be used with more 
confidence. At first one is tempted to say that the virial ra- 
dius should provide a positive answer, since at this radius 
non-gravitational processes are thought to be less important. 
However, there are two aspects that do not allow us to make 
definite conclusions on this point: a) the observations of high- 
redshift clusters do not extend up to the virial radius and so the 
gas masses at this region were obtained by extrapolating the 
y8-model fitted on the observed parts of the clusters, and b) the 
maximum likelihood values for Q,„ at this radius are too high 
(Q„, ~ 0.7), which translates to a significant deviation from 
the concordance model supported by other observational data 
such as CMB, Type la supernovae, and baryonic oscillations 
(Spergel et al. 20 03] lAstier et al. 20061 lEisenstein et al. 20051 
[Blanchard et al. 20061 1. Observations of X-ray emissivity in 
distant clusters up to the virial radius will be critical to pro- 
viding definitive answers on the first aspect. Despite this, the 
use of gas mass fraction evolution (or lack of it) in clusters as 
a precision method for constraining cosmological parameters 
seems to be seriously compromised. 
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